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Abstract 


The SARS-associated human coronavirus (SARS-HCoV) is a newly described, emerging virus conclusively established as the etiologic 
agent of the severe acute respiratory syndrome (SARS). This study presents a single-tube RT-PCR assay that can detect with high analytical 
sensitivity the SARS-HCoV, as well as several other coronaviruses including other known human respiratory coronaviruses (HCoV-OC43 and 
HCoV-229E). Species identification is provided by sequencing the amplicon, although a rapid screening test by restriction enzyme analysis 
has proved to be very useful for the analysis of samples obtained during the SARS outbreak in Toronto, Canada. 


© 2004 Elsevier B.V. All rights reserved. 
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1. Introduction 


Severe acute respiratory syndrome (SARS) is a newly de- 
scribed, emerging infectious disease presenting with high 
fever and atypical pneumonia. The disease first appeared in 
the Guangdong province of China and quickly spread to other 
countries; outbreaks were reported in mainland China, Viet- 
nam, Singapore, Hong Kong, Taiwan and Canada (Drazen, 
2003; Lee et al., 2003; Poutanen et al., 2003; Tsang et al., 
2003; Vu et al., 2003; WHO, 2003a). A new, previously 
unknown coronavirus was first isolated from patients with 
SARS at the University of Hong-Kong (Peiris et al., 2003a; 
WHO, 2003a); this was quickly followed by reports from 
several other laboratories of the demonstration of a new coro- 
navirus in samples from patients with SARS (Drosten et al., 
2003; Ksiazek et al., 2003; Poutanen et al., 2003; WHO, 
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2003a). The complete genomic sequence of the virus has 
been determined (Marra et al., 2003; Rota et al., 2003) and 
confirmed the classification of the SARS-associated agent 
as a new coronavirus, distinct from the previously known 
groups of coronavirus (Marra et al., 2003; Rota et al., 2003) 
but with some relationship with group 2 (Snijder et al., 2003). 
Experimental infection of cynomolgous monkeys permitted 
the fulfillment of Koch’s postulates for this pathogen and 
firmly established the SARS-associated human coronavirus 
(SARS-HCoV) as the etiologic agent of SARS (Fouchier 
et al., 2003). 

This study presents a RT-PCR protocol that allows for 
the detection of the SARS-HCoV in clinical samples. This 
single-tube RT-PCR is based on consensus primers target- 
ing conserved regions of coronavirus genome sequences and 
allows for the detection and species identification of sev- 
eral coronaviruses including SARS-HCoV, with high an- 
alytical sensitivity. This assay is expected to be helpful 
in fulfilling the World Health Organization laboratory case 
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definition of SARS, which recommends the use of two dif- 
ferent PCR assays (WHO, 2003b). In addition, even if the 
SARS-HCoV does not reappear in the human population, 
this assay will still be useful, since it provides a diagnos- 
tic test for the respiratory coronaviruses HCoV-229E and 
HCoV-0C43. 


2. Material and methods 
2.1. SARS coronavirus RNA 


RNA was extracted with TRIzol reagent (Life Technolo- 
gies), as per manufacturer’s recommendations, from a lung 
biopsy from a patient with SARS. The RNA pellet was re- 
suspended in 30 wl of 10 mM of dithiotreitol with 5% (v/v) 
of RNasin (20-40 U/l, Promega), serially diluted, aliquoted 
and frozen at —80 °C. Quantification of the RNA serial dilu- 
tion was achieved by measuring the amount of SARS-HCoV 
RNA in the aliquots with the RealArt HPA coronavirus RT- 
PCR (Artus GmbH, Hamburg, Germany). 


2.2. RNA from human coronaviruses OC43 and 229E 


Human coronaviruses OC43 and 229E strains were ini- 
tially obtained from ATCC and passaged in cell culture as 
described (Sizun et al., 1998). Stocks of viruses were titrated 
as described (Sizun et al., 1998). RNA was extracted from 
the titrated stocks using TRIzol (Life Technologies) as per 
the manufacturer’s recommendations and serially diluted. 


2.3. RNA from infectious bronchitis virus 


Tissue culture-adapted infectious bronchitis virus (IBV) 
stocks, strain Massachussets 41 and strain Baudette, were 
obtained from Dr. E. Nagy, Ontario Veterinary College, Uni- 
versity of Guelph, Ontario. RNA was extracted and serially 
diluted as for HCoV-OC43 and HCoV-229E. 


2.4. Clinical samples 


During the SARS outbreak in Toronto (2003), several sam- 
ples from patients with probable or suspected SARS were 
referred to our laboratory for SARS-HCoV detection by RT- 
PCR. Of these, 44 were assayed using the new primers de- 
scribed here, including 7 lung biopsies, 2 bronchoalveolar 
lavages, 7 nasopharyngeal swabs, 3 throat swabs, 2 pleu- 
ral fluid samples, 2 endotracheal aspirates, 4 blood samples, 
11 stool samples and 6 urine samples. RNA was extracted 
from lung tissue using the RNeasy kit (Qiagen) as per 
the manufacturer’s recommendations. For all other samples, 
RNA was extracted using the guanidium thiocyanate buffer 
(GTC) extraction methods, as described (Johnson et al., 
2000). 


2.5. Primers 


Primers were designed to target segments of the pol Ib 
coding region conserved among several species of coron- 
avirus (Fig. 1). We used the sense primer CORO1 (5’-TGA 
TGG GTT GGG ACT ATC CTA AAT GTG A-3’) and the an- 
tisense primer CORO2 (5’-GTA GTT GCA TCA CCG GAA 
GTT GTGCCA CC-3’), homologous to the segments [15211, 
15238] and [15402, 15430] of the sequence of the Tor2 strain 
of SARS-HCoV (GenBank Accession number AY274119), 
respectively. 


2.6. RT-PCR 


This was done using the Qiagen one-step RT-PCR kit. 
Each reaction was performed in a 0.6 ml tube (Diamed PRE 
050) in a total volume of 50 wl overlaid with 50 wl of mineral 
oil. Each reaction mix contains 10 wl of 5x Qiagen buffer, 
2 wl of dNTP mix (each dNTP at 10 mM concentration), 3 pl 
of each primer (10 1M stock), 20 wl of molecular grade dou- 
ble distilled water (ddH20O) and 2 wl of the one-step enzyme 
mix (Qiagen). The master mix was then aliquoted in tubes, to 
which 10 wl of template RNA was added. The PCR thermal 
cycling was done on a Stratagene Robocycler 40, with an ini- 
tial incubation at 50 °C for 30 min, followed by an incubation 
at 95 °C for 15 min, and 40 cycles consisting of denaturation 
at 94°C for | min, annealing at 56°C for | min, and elonga- 
tion at 72°C for | min. Extensive precautions against PCR 
contamination, as previously described (Johnson et al., 2000), 
were strictly observed. Positive and negative controls, as well 
as extraction controls and controls for PCR inhibition, were 
set up essentially as described (Johnson et al, 2000). 


2.7. Electrophoresis analysis of the amplicons 


A 10pl volume of each reaction was submitted to 
electrophoresis on 2% agarose or 4% NuSieve 3:1 (Biowhit- 
taker Molecular Applications, Rockland, USA) gels contain- 
ing ethidium bromide. The gels were visualized on a UV 
transilluminator and photographed. 


2.8. Restriction enzyme analysis 


Each reaction mixture in which the expected 220 bp am- 
plicon was detected was digested with the restriction enzyme 
Alul, which cuts within the sequence AGCT. The digestion 
reaction consisted of 10 wl of the RT-PCR mixture, 1.5 jl of 
the 10x enzyme buffer, 1 wl of AluI (10 U/ul, New England 
Biolabs) and 2.5 wl of ddH2O for a total of 15 wl. The reac- 
tion was incubated at 37°C for 1h and analyzed by agarose 
gel electrophoresis. 


2.9. Sequencing 


Amplicons were submitted to automated sequencing, for 
both strands, using the PCR primers as sequencing primers. 
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Fig. 1. Alignment calculated from the nucleotide sequence of segments homologous to the segment [15211, 15430] of SARS-HCoV strain TOR2 (GenBank 
Accession number AY274119) from 12 other species of coronaviruses. HCoV-NL63 (AY567487): human coronavirus NL63. HCoV-229E (NV_002645): human 
respiratory coronavirus 229E. PEDV (NC_003436): porcine epidemic diarrhea virus. FIPV (AF124987): feline infectious peritonitis virus. CCoV (AF124986): 
canine coronavirus. TGEV (AJ271965): porcine transmissible gastroenteritis virus. SDAV (AF124990): sialodacryoadenitis virus. MHV (NC_001846): murine 
hepatitis virus. BCoV (NC_003045): bovine coronavirus. HCoV-OC43 (AF124989): human respiratory coronavirus OC43. IBV (NC_001451): avian infectious 
bronchitis virus. TCoV (AF124991): turkey coronavirus. Also aligned with the coronaviruses sequences are the sequences of the sense primer (CORO1) and 
the complement of the antisense primer (CORO2rc). Dots denote homology with the consensus sequence on top. The sequence alignment was calculated using 
ClustalX for Windows version 1.81 (Thompson et al., 1997) and edited with Genedoc version 2.3 for Windows (Nicholas K.B., 1997). 


Sequencing was performed by the DNA Sequencing Fa- 
cility, Centre for Applied Genomics, Hospital for Sick 
Children. 


2.10. Sequences analysis 


Sequence editing and analysis were done using the pro- 
grams Generunner for Windows version 3.05 (Hasting soft- 
ware). Sequence alignments were calculated using ClustalX 
for Windows version 1.81 (Thompson et al., 1997), with the 
default parameters for gap opening and gap extension. Phy- 
logenetic trees were inferred using Treecon for Windows ver- 
sion 1.3.b (Van de Peer and De Wachter, 1994), using a dis- 
tance method; in brief, the distance was calculated without 


correction; the tree topology was inferred with the neighbor 
joining method and the trees re-rooted at the internode. Boot- 
strap analyses were done with 1000 replicates. 


3. Results 
3.1. Sequence alignments and primer design 


Initially, an amino acid sequence alignment of the 
pol 1b open-reading frame (ORF) was used to identify 
well-conserved regions across 12 species of coronaviruses 
comprising the SARS-HCoV and members of the three previ- 
ously described groups of coronaviruses. The corresponding 
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Fig. 2. (A) Phylogenetic tree calculated from the nucleotide sequence of segments homologous to the segment [15211, 15430] of SARS-HCoV strain TOR2 
from the 13 species of coronaviruses in Fig. 1. 1. The boostrap values are displayed (as percentages) above the node. (B) Phylogenetic tree calculated from the 
amino acid sequences corresponding to the nucleotide sequences in (A). The boostrap values are displayed (as percentages) above the node. The phylogenetic 
trees were drawn using the program Treecon for Windows version 1.3.b (Van de Peer and De Wachter, 1994). 


nucleotides alignment was then used to design PCR primers 
targeting these well-conserved regions. Fig. | illustrates the 
alignment of the 220 nt segments of these 12 coronavirus se- 
quences homologous to the segment [15211, 15430] of the 
Tor2 strain of SARS-HCoV, as well as the corresponding 
sequence of the newly described human coronavirus HCoV- 
NL63 (Fouchier et al., 2004; Van der Hoek et al., 2004), along 
with the sequence of the primers. 

Fig. 2 illustrates the phylogenetic trees of the nucleotide 
sequences (Fig. 1A) and amino acid sequences (Fig. 2B) in 
that segment, demonstrating that there are enough sequence 
diversity in the region bracketed by the primers to allow for 
unambiguous species identification. 


3.2. Detection of coronavirus RNA 


Using as templates the RNA of SARS-HCoV, HCoV-229E 
and HCoV-OC43, as well as RNA from two strains of IBV, 
we could readily and reproducibly obtain the expected 220 bp 
amplicon. 


3.3. Sensitivity and specificity 


Aliquots of a 10-fold serial RNA dilution prepared from 
a lung biopsy sample of a patient with SARS (see Section 2) 
were used to compare our assay with the RealArt HPA coro- 
navirus RT-PCR (Artus GmbH). The results are displayed 
in Table 1 and show that our assay has essentially the same 
sensitivity as that of the RealArt HPA coronavirus RT-PCR 
assay; it is estimated that the analytical sensitivity is between 
1 and 10 genome copies. 

For HCoV-OC43, viral RNA could be detected from as 
little as 6 x 1072 TCIDs (tissue culture infective dose 50); 
for HCoV-229E, viral RNA could be detected from as little as 
6 x 10-3 TCIDs0. For all coronaviruses tested in this study 
the species was verified by sequencing; in all cases the se- 
quence internal to the primers were identical to the expected 


sequences. Based on the expected sequences of the ampli- 
cons a rapid screening test for species identification for the 
human coronaviruses was designed using the A/uI restriction 
enzyme. As illustrated on Fig. 3, the SARS-HCoV ampli- 
con is cut into 126-bp and 94-bp fragments, the HCoV-229E 
amplicon into 101-bp, 86-bp and 33-bp fragments, and the 
HCoV-OC43 amplicon is not cut by Alul; based on the refer- 
ence sequence, the amplicon of IBV would not contain a Alul 
site. Nucleic acids from human tissue, or from other respira- 
tory viruses whose presence was detected in various samples 
(including influenza A, respiratory syncytial virus, parain- 
fluenza 3 and adenoviruses), did not generate amplicons of 
220 bp. 


3.4. Clinical samples 


From the clinical samples listed above, the presence of 
the SARS-HCoV RNA was demonstrated in one lung biopsy, 
two bronchoalveolar lavages, two nasopharyngeal swabs, two 
urine samples and one throat swab. 


Table 1 
Comparison between our RT-PCR assay and the RealArt HPA coronavirus 
RT-PCR (Artus GmbH) 


Dilution RT-PCR with CORO1 


RealArt HPA coronavirus RT-PCR 


and CORO2 (genome copies) 
10-3 + 33000 
10-4 + 2480 
10-5 + 89 
10-6 + 19 
10-7 + 3 
10-8 = - 


Column two shows the result (scored as positive or negative) obtained with 
our assay using 10 pl of the indicated RNA dilution (see Section 2). For 
each dilution the assay was done on two different aliquots and the results 
were in agreement. Column three shows the average (on two experiments) 
of the measured number of genome copies in 5 pl of the RNA dilution with 
RealArt HPA coronavirus RT-PCR. 


D. Adachi et al. / Journal of Virological Methods 122 (2004) 29-36 33 


SARS — HCoV 


= 
= 


MW 


HCoV — 229E 
HCoV — OC43 


MW 


1 2 3 


4 5 6 


Fig. 3. Alul digestion of amplicons from SARS-HCoV, HCoV-229E and HCoV-OC43. Lanes 1, 3 and 5 (MW): molecular weight markers (Amplisize ladder 
50-2000 bp, Bio-Rad Laboratories, Mississauga, Canada); Lane 2: SARS-HCoV (predicted fragments: 126 bp, 94 bp); Lane 4: HCoV-229E (predicted fragments: 


101 bp, 86 bp, 33 bp); Lane 6: HCoV-OC43 (predicted fragment: 220 bp). 
4. Discussion 


This study presents a single-tube RT-PCR assay designed 
to detect several species of coronaviruses, including the 
SARS coronaviruses, by using consensus primers designed 
from an alignment of the sequences of a region in the pol 1b 
ORF of 12 species of coronaviruses. A previous version of 
these primers, based on an alignment of four species of coron- 
aviruses and of the Berne torovirus, was used to demonstrate 
the presence of a new coronavirus in patients with SARS 
in Canada (Kumar et al., 2003; Poutanen et al., 2003) at a 
time when the complete sequence of the SARS coronavirus 
was still unknown. The use of the revised primers and the 
optimization of the RT-PCR has led to an improvement in 
analytical sensitivity of three orders of magnitude compared 
to the previous assay. Its sensitivity has proven sufficient for 
the detection of the SARS-HCoV in several clinical samples. 

The assay was designed to be broadly reactive with the 
genome of many coronavirus species; it is demonstrated 
here that it can detect coronaviruses from all four known 


groups of coronaviruses, including the HCoV-229E and 
HCoV-O0C43, which are well-recognized human pathogens 
(Holmes, 2001). Although not all the known coronavirus 
species were tested, it is expected that the primer pair can 
amplify with high analytical sensitivity the targeted genome 
segment from the 12 species used in the design of the primers, 
with the possible exception of the porcine epidemic diarrhea 
virus (PEDV) whose sequence has a mismatch at the third 
nt from the 3’ end of the primer CORO1 (Fig. 1). The cur- 
rent assay can detect with high sensitivity three coronaviruses 
known to infect humans, and sequencing of the amplicons 
would unambiguously identify the species (Fig. 2). Although 
there has not been yet an opportunity to test this assay on 
isolates of the newly described HCoV-NL63 (Van der Hoek 
et al., 2004), it is predicted, based on the sequence alignment 
in Fig. 1, that our RT-PCR assay would also detect HCoV- 
NL63, and sequencing of the amplicon followed by phyloge- 
netic analysis (Fig. 2) would readily confirm the identity of 
the virus. Based on the published sequence of HCoV-NL63, 
digestion of the amplicon by AluI would yield a distinct 
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pattern consisting of bands of 187 bp and 33 bp. Since the 
sequences targeted by the primers CORO! and CORO2 are 
conserved across several coronavirus species (and this is par- 
ticularly clear at the amino acid level), this suggests strong 
selective pressure on that part of the genome. In turn, this 
suggests that the assay would still be able to detect mutant 
or variant strains of SARS-HCoV, as mutations would likely 
not occur in the regions targeted by the primers. To date, the 
various isolates of SARS-HCoV have all displayed a great 
genetic homogeneity (Ruan et al., 2003), which argues for 
the emergence of SARS-HCoV in the human population only 
very recently. However, should the virus emerge again and be 
allowed to propagate, its genetic diversity may well increase 
(although the genomes of many known coronaviruses appear 
relatively stable). Of note, a BLAST search in the GenBank 
database using the 220-nt segment of SARS-HCoV (Fig. 1) 
showed a 100% homology with 105 sequences of SARS- 
HCoV deposited to date, including representative of all clus- 
ters (Guan et al., 2004), and an homology greater than 99% 
(219/220) with two other sequences of SARS-HCoV. As well, 
the recent SARS outbreak suggests that other coronaviruses 
may emerge as new human pathogens. Again, assuming that 
the conservation of the sequences targeted by the primers 
stems from strong evolutionary constraints, the RT-PCR as- 
say described here would very likely detect such an emerging 
coronavirus. 

In this study a simple restriction enzyme analysis as a 
screening test was shown to be very helpful in identifying 
the SARS-HCoV. Indeed, in all the samples that tested pos- 
itive for SARS-HCoV in our laboratory during the Toronto 
outbreak, the expected restriction pattern was demonstrated 
(and based on sequences deposited in GenBank this pattern 
would have been demonstrated in all isolates of SARS-HCoV 
sequenced to date). However, because RNA virus genomes 
are much more variable than DNA virus genomes, one should 
not rely completely on such an analysis and the failure to ob- 
tain the expected digestion fragments (or if the amplicon is 
not cut, as is the case for HCoV-OC43) should prompt the 
determination of the sequence of the amplicon. As an illus- 
trative case in point, whereas the amplicon from the IBV 
strain Baudette is not cleaved by AluI, the amplicon from 
strain Massachussetts 41 is in fact cleaved by Alul, as verified 
by both sequencing and digestion of the amplicon (data not 
shown). After sequencing the amplicon, a comparison of the 
sequence to that of several coronaviruses in a phylogenetic 
analysis, as in Fig. 2, should permit a definitive identification. 
Indeed, should a new coronavirus emerge, it appears likely 
that it could be detected and identified as a new virus by using 
our RT-PCR assay followed by sequencing and phylogenetic 
analysis. 

Since the end of the global SARS outbreak in 2003, a to- 
tal of three laboratory confirmed sporadic cases of SARS in 
Guangdong have occurred (WHO, 2004), which however did 
not spread further, as well as three laboratory acquired inci- 
dents (Normile, 2004). While surveillance for SARS is still 
recommended (WHO, 2003b) there is great concern about 


the negative impact of a false positive laboratory test, and the 
World Health Organization has recommended protocols for 
the laboratory diagnosis of SARS in the post outbreak period 
(WHO, 2003b), which includes the use of two different RT- 
PCR assays. The assay described here provides an alternative 
test that can be helpful in that regard. Furthermore, because 
our assay detects several coronaviruses, it is quite feasible 
to use RNA from a coronavirus other than SARS-HCoV as 
a positive control. In fact, our laboratory now uses as posi- 
tive control RNA transcribed from the cloned amplicon of the 
IBV; IBV has never been implicated in human infections. The 
use of IBV as a positive control removes a potential source 
of PCR contamination by SARS-HCoV genetic material in 
the laboratory. 

In this study, the presence of the SARS-HCoV was demon- 
strated in a variety of clinical samples, which illustrates the 
usefulness of the assay. The number of positive samples re- 
ported here is too small to draw firm conclusions but taking 
into account the distinction between probable and suspected 
cases of SARS, as well as previous observations made ear- 
lier in our laboratory during the outbreak, it appears that the 
SARS-HCoV is more likely to be demonstrated in lower res- 
piratory tract samples, such as broncho-alveolar lavages or 
lung biopsies rather than in a nasopharyngeal swab; this par- 
allels the findings of other reports (Peiris et al., 2003b; Tang 
et al., 2004) and it has been suggested that the SARS-HCoV 
appears at a high titer in the nasopharynx relatively late in 
the disease (Peiris et al., 2003b). Thus, the negative RT-PCR 
results on such samples do not reflect a poor analytical sensi- 
tivity of the assay but rather the biology of the virus. A system- 
atic study of various samples obtained throughout the course 
of the disease in several patients will be required to determine 
the best testing algorithm. Certainly in the control of the out- 
break in Toronto, carefully established epidemiological links 
between patients fitting the case definition seemed the most 
effective method, but laboratory confirmation of some cases 
in the chain also proved to be very helpful (Wallington et al., 
2003). 

In summary, a single-tube RT-PCR assay for the detec- 
tion of human coronaviruses, including the SARS-HCoV, has 
been developed. It is expected to remain effective with pos- 
sible SARS-HCoV mutants and to be of value should other 
new coronaviruses emerge as human pathogens. 
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